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Abstract 

In this short review it is suggested that the relationship between 
third- and fourth-order moments of turbulence in the atmospheric 
boundary layer depends on stability. This can explain some differences 
among datasets, and provides a key point for modelling improvement. 



1 Introduction 



Data on the third- and fourth-order moments of turbulent velocities in bound- 
ary layers have been collected for many years, not only in the at mospheric 



boundary layer (ABL), where their evaluation is somewhat difficult ()Lenschow et al 
1994^ . but also in laboratory experiments, where the relations hip between 



odd-o rder moments and their next even-order was studied by iDurst et al. 
fll987ri. for instance. 



Cheng et all. Isubmittedh 



Some recent papers (jFerrero and R,accaL I2C 
discuss the role of high order moments in parameterisations of the ABL, 
often referred to as non-local models in that space derivatives of high order 
moments are used to model lower orders. 
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The role of third order moments has long been recognised, at least in 
the Convective Boundary Layer (CBL), as connected to non-local trans- 
port properties ( Wvngaard and Weill . 1991 ). As regards the fourth-order 
mo ments, in the lack of any f urther information, the Millionchikov hypothe- 
sis fjMonin and YaglomLll97lL p. 241) is usually invoked. This means that for 
vanishing skewness that tends to zero, the normalised fourth order moments 
(kurtosis) tend to their Gaussian values. 

The question is whether the available data support this hypothesis, to 
what extent and under which conditions. 



2 Review of data and parameterisations 



In the study of high-order moments, particular interest is devoted to the nor 
mahsed mor nents, i.e., skewness (.S ) and kurtosis {K) and r elationships be- 
tween them. Maurizi and Tampieril ()l999l ) , and subsequently Tampieri et al. 



( 2000l ). collected and analysed data from literature on horizontal and vertical 
velocity skewness and kurtosis, considering a variety of turbulence genera- 
tion mechanisms. Furthermore, recent papers foc us attention in particu- 
lar on the CBL: new measurements from aircraft JOrvanik and H^rtmani 
2OO2I GH herea f ter) a nd with remote sensing acoustic techniques (SODAR) 
(|Alberghi eraP . 120021 AMT hereafter) have been presented 

D ifferent simple parameterisations were propose d for the S-K relation- 
ship. iTampieri et al . (2000), first, and subsequentlv iMaurizi and Lorenzani 
( 20011 ). proposed a parameterisatio n based on the observation that a statis- 
tical limit exists in the S-K space flKendall and Stuard . Il977l l. namely 



K> K, 



lim 



5^ + 1 



(1) 



This limit shapes the structure of the S-K. Thus, the pair (S', K) with 
K = K{S^ + 1)~^ can be taken as the natural coordinate system for the 
S-K space. The simplest model (zero order) based on this observation can 
be built assuming a constant K and, therefore, 

K = ao{S^ + 1) . (2) 

Fitting Equation Q to data, Tampieri et al.l ( 2000l ) found, for the ver- 
tical velocity component, = 3.3 for shear dominated boundary layers 
(see Figure HI) and a n = 2.5 in the CBL (AMT), adding new data to the 
Tampieri et al. (|2000>) dataset, confirmed the CBL result, giving ao = 2.4. 
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Figure 1: Skewness and kurtosis for shear dominated boundary layers 
( Tampieri et al. . 2000[ ). Continuous line represents the statistical limit, while 
the dashed line is Equation Q with = 3.3. (reprinted courtesy of ... ) 
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GH found that a mass-flux assumption for the CBL exactly results in 
K = Kiijn (see Equation (jT))) as a relationship between S and K. In fact, 
it is known that this relationship only holds for two-value processes. They 
used, generalisation, the form 

K = aoif3S^ + 1) . (3) 



It is worth noting that Lewis et'aD (jl997l ) proposed a relatio nship of the same 



form f or concentration data, based on the ideas expressed bv lChatwin and Sullivan 



[t can be observed that, for \S\ < 1, Equation Q is consistent with 



a second order model in the (5*, K) space, namely 

K = ao + aiS + 028'^ , (4) 
with ai = and 0:2 = aoiP — !)■ Because the d a taset of GH dataset shows 



strictly leptokurtic cases (see also ICheng et all Isubmittedf ) . they selected 



ao = 3, ai = and 0:2 = —2 in Equation (jH), and thus assumed Gaussianity 
for symmetric distributions. 

This assumption is not in agreement with the AMT dataset, which shows 
the presence of a great deal of data below K = 3. However the data in GH 
refer to a particular case (cold air outbreak over the ocean) and fall within 
the area covered by the AMT dataset, which collects measurements from 
presumably different environmental conditions. 

The question remains as to whether this subset is simply incomplete with 
respect to a full convective behaviour, or whether it reflects the fact that 
the production mechanism determines the S-K relationship. In fact, while 
the presented data are nominally taken in convective conditions, a concur- 
rence of different production mechanisms (shear and buoyancy) is actually 
expected. For instance, the Monin-Obukhov similarity theory states that 
for z -C |L| shear dominates over buoyancy, and it is presumable that data 
taken at different z/L present a different balance between the two mecha- 
nisms. Furthermore, the properties of measurements taken across the CBL 
depth as a whole could depend on \L\/zi, which gives a measure of the frac- 
tion of CBL where buoyancy cannot be considered as the only production 
mechanism. 

The two values ao = 2.4 to oq = 3.3 in Equation ^ suggest that the 
Gaussian case {K ^ 3 as 5 ^ 0) may occur for the vertical turbulent velocity 
as a transition between the shear and convective production mechanisms (for 
instance, in the CBL at 2; < Lmo)- 
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Figure 2: Data of S and K for vertical velocity measured in turbulent bound- 
ary layer over an obstacle at three different positions x with respect to the 
obstacle top. Full triangle: x = (Ric > 0), open circle: x = L (Ric = 0), 
open triangle: x = — oo (Ric < 0). 

This observation can be expressed by letting be a function of, at least, 
the Richardson number Ri, thus re-writing Equation Q as 

K = aoiRi)iS^ + I) (5) 



with the constraint ao(Ri) > 1. It can be argued that a n increases as z/Lmq 
increases from negative values to zero, consistently with lAnfossi et al.l (|l997 . 
their Table 1). 



Rc<0 


Rc = 


Rc>0 


2.8 


3.0 


3.2 



Table 1: Minimum value of K measured in three different stability conditions. 

3 More on stability effects 

It would be interesting to investigate this issue further if certain details on 
the measurements were available. In order to overcome the unavailability of 
such information, we consider a flow over a simple obstacle, and use the anal- 
ogy between buoyancy and streamline curvature (see, e.g., Bradshaw . 19691 : 



Baskaran et al.L 119911 1 to investigate on the role of stability in determining 
the character of the S-K relationship. 

In a neutral turbulent flow with streamline curvature, such as a flow over 



a hill , it is possible to define a curvature Richardson number Ric ( Wvngaardl . 



19671) as 

2U 

= du^ u 

dz ^ R 

where U is the mean velocity module and R is the streamline curvature ra- 
dius. Positive Ric corresponds to a dumping term in the turbulent energy 
budget (e.g., on the hill top) expressed in streamhne coordinates, while neg- 
ative Ric represents an unstable contribution, e.g., near the hill base. 

An experiment carried out in the Enfio "A" wind tunnel focused attention 
on third- and fourth-order turbulence statistics . A turbulent boundary layer 
was generated upstream of a sinusoidal, two-dimensional hill with aspect 
ratio H/L = 5, where H is the obstacle height and L is half the total width. 
Turbulence was measured with a hot-wire anemometer at 7 different positions 
streamwise, at 20 vertical levels. Measurements were recorded for times long 
enough to produce reliable statistics up to the fourth order. 

Figure 121 shows data of 5* and K for three profiles measured far upstream, 
on the hill top and on the lee side at the hill foot. Those regions correspond 
to Ric = 0, Ric > and Ric < 0, respectively. Apart from the specific 
structure revealed by the measurements in the unstable case for S in the 
range 0.2-0.4, the form of the S-K relationship is similar for the three cases 
and exhibits a minimum value. 
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We use this minimum as a rough indicator of the different features of 
the S-K relationship. This seems a reasonable parameter in that it could 
discriminate data reported by GH from other data in the AMT dataset. 

The minimum measured values of K are report ed in Table H] for th ree 
different cases of Ric estimated from a similar flow ( Maurizi et al.l . ll997 ). 

Although measured profiles are largely influenced by the specific dynamics 
of flow over the obstacle and, in particular, by non local equilibrium in the 
lee side, results are in qualitative agreement with the proposed dependence 
of expansion coefficients on stability (Equation (0)). 



4 Conclusions 

The S-K relationship is modelled by an expansion in an appropriate space 
and some parameterisations in literature have been reduced to this scheme. 

In this frame, comparing data for shear- and convective-dominated bound- 
ary layers, it is found that model constants should depend on stability. In 
particular, the Gaussian case can occur in intermediate situations between 
the purely shear- and purely convective-dominated boundary layer. 

Furthermore, an examination of the differences between two datasets for 
the convective case suggests that that there could be a variety of behaviours 
for different balances between shear and buoyancy production. 

As an example, some data measured in a wind tunnel experiment have 
shown the possibility that the suggested dependence can be confirmed. It is 
worth pointing out that this dataset if far from exhaustive with respect to 
the problem studied and it has been used merely as a indicator. 

Further measurements of high-order moments of wind velocity in really 
different stability conditions are required in order to provide a more quanti- 
tative response to the problem. 
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